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Taking  into  account  inner  and  external  multi-irreversibilities,  a  complete  numerical  model  of  commercial 
thermoelectric  generator  with  finned  heat  exchangers  is  established  by  combining  thermodynamics  with 
heat  transfer  theory.  A  significant  novelty  is  that  physical  properties,  geometric  dimensions,  temperature 
parameters  and  flow  parameters  are  all  considered  in  the  model.  The  inner  effects  include  Seebeck  effect, 
Fourier  effect.  Joule  effect  and  Thomson  effect.  The  irreversibilities  include  the  heat  transfer  through  the 
air  gap  (proposed  and  evaluated  first  time),  the  thermal  and  electrical  resistance  of  the  conducting  strips, 
and  the  multiform  external  thermal  resistances.  Based  on  the  numerical  model,  the  performances  of 
a  typical  commercial  thermoelectric  generator  are  simulated.  Hot  water  at  60—100  °C  and  cold  water  at 
27  °C  are  employed  as  heat  source  and  sink  of  the  generator  module  which  consists  of  127  thermoelectric 
elements.  The  results  show  that  the  maximum  power  output  of  0.13  W  and  the  maximum  efficiency  of 
0.87%  are  available  from  the  generator.  The  open  circuit  voltage  is  1.80  V  and  the  short  circuit  current  is 
0.28  A,  respectively.  The  effects  of  external  irreversibilities  on  the  performance  of  the  thermoelectric 
generator  are  analyzed  by  comparing  this  irreversible  model  with  the  exo-reversible  model.  The 
numerical  model  and  calculation  method  can  be  applied  to  the  performance  prediction  and  optimization 
of  thermoelectric  generators  with  finned  heat  exchangers.  The  simulation  results  can  be  used  as  feasi¬ 
bility  and  effectiveness  reference  by  employing  low-grade  energy  or  waste  heat  for  power  generation. 

©  2011  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

There  are  several  processes  that  can  create  electrical  current 
from  a  thermal  gradient:  thermoelectric  [1],  thermionic,  thermo¬ 
magnetism,  ferroelectricity,  and  the  Nernst  effect.  Thermoelectric 
conversion  is  the  most  effective  of  these  processes  [2,3].  Moreover, 
the  most  successful  approach  employing  indirect  conversion  of 
alpha  or  beta  emissions  into  electricity  is  using  thermoelectric 
effect.  That  is  called  a  radioisotopic  thermoelectric  generator  (RTG) 
[4]  which  is  used  in  Micro-electromechanical  Systems  (MEMS) 
frequently.  As  a  low-carbon  energy  technology,  thermoelectrics 
have  been  used  in  military,  aerospace,  instrument,  and  industrial  or 
commercial  products,  as  a  power  generation  for  specific  purposes. 
More  widespread  use  of  thermoelectric  requires  not  only  improving 
the  intrinsic  energy-conversion  efficiency  of  the  materials  [5,6],  but 
also  modeling  and  analyzing  of  thermoelectric  device  system  [7,8], 

In  general,  conventional  non-equilibrium  thermodynamics 
[9—11  ]  is  used  to  analyze  the  performance  of  one  or  multi-element 
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thermoelectric  generators.  Considering  the  inner  structure  of 
a  thermoelectric  generator,  a  significant  increase  in  the  electrical 
power  output  from  a  module  can  be  achieved  by  modifying  the 
geometry  of  the  thermoelectric  elements  [12,13].  Rowe  investi¬ 
gated  the  efficiency  of  a  single  couple  solar  powered  thermoelectric 
generator  [14]  and  reviewed  US  applications  of  nuclear-powered 
thermoelectric  generators  in  space  [15],  Sisman  and  Yavuz  [16] 
analyzed  the  effect  of  Joule  losses  on  the  total  efficiency  of  a  ther¬ 
moelectric  power-cycle.  Chen  et  al.  [17]  investigated  the  influence 
of  Thomson  effect  on  the  maximum  power  output  and  maximum 
efficiency  of  a  thermoelectric  generator.  Xuan  et  al.  [18]  employed 
a  phenomenological  model  to  study  the  effects  of  internal  and/or 
external  interface  layers  on  thermoelectric  devices  performance. 

Thermoelectric  devices  cannot  be  used  independently.  They 
should  be  connected  with  heat  exchangers  to  absorb  and  dissipate 
heat  [19-21],  The  theory  of  finite  time  thermodynamics  (FTT)  or 
entropy-generation  minimization  [22-31  ]  is  a  powerful  tool  for  the 
performance  analysis  and  optimization  of  practical  thermodynamic 
processes  and  devices  with  heat  transfer.  Some  authors  have 
investigated  the  performances  of  thermoelectric  generators  using 
a  combination  of  finite  time  thermodynamics  and  non-equilibrium 
thermodynamics.  Considering  the  external  irreversibilities  of 
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Nomenclature 

v 

kinematic  viscosity  (m2  s  1) 

2 

(O 

intermediate  variable 

Aa 

thermal  diffusivity  (m2/s) 

Subscripts 

c 

heat  capacity  ratio  0  kg  1  K-1) 

a 

side 

ej 

electrical  current  density  (A  m  2) 

b 

base,  inner  surface 

F 

intermediate  variable 

c 

cold  junction,  channel 

g 

gravitational  acceleration  (m  s-2) 

cd 

heat  conduction 

H 

height  (m) 

cp 

ceramic  plate 

h 

heat  transfer  coefficient  (W  m  2  K-1) 

cu 

copper  conducting  strip 

/ 

electrical  current  output  (A) 

cv 

heat  convection 

K 

thermal  conductance  (W  K-1) 

e 

external 

k 

thermal  conductivity  (W  m-1  K-1) 

ex 

exchanger 

L 

length  (m) 

er 

external  reversible 

m 

intermediate  variable 

f 

fin,  fluid 

N 

number 

g 

generation,  air  gap 

P 

power  output  (W) 

H 

heat  source 

P 

perimeter  (m) 

h 

hot  junction 

Q 

heat  flow  rate  (W) 

i 

exterior  surface 

R 

electrical  resistance  (Q),  thermal  resistance  (K  W-1) 

in 

input 

S 

entropy  generation  rate  (W  K-1) 

J 

Joule  heat 

T 

temperature  (K) 

K 

heat  conduction  leakage 

U 

voltage  output  (V) 

n 

N-type  semiconductor  leg 

V 

velocity  (m  s-1) 

out 

output 

x 

coordinate 

P 

P-type  semiconductor  leg 

Z 

figure  of  merit  of  a  thermoelectric  element  Z  =  a 2I(KR) 

rd 

heat  radiation 

T 

total 

Creek  letters 

te 

thermoelectric  element 

a 

Seebeck  coefficient  (V  K_1) 

0 

coefficient  of  cubical  expansion  (K-1) 

Superscripts 

A 

difference 

P 

at  maximum  power  output 

8 

thickness  (m) 

V 

at  maximum  efficiency 

e 

blackness 

cp 

kinetic  viscosity  (kg  m  1  s~ ') 

Abbreviation 

V 

Efficiency 

Gr 

Grashof  number  Gr  =  g^\TL3lv2 

fi 

Thomson  coefficient  (V  K-2) 

Nu 

Nusselt  number  Nu  =  hH/kf 

e 

packed  density 

Pr 

Prandtl  number  Pr  =  v/a 

p 

electrical  resistivity  (Qm) 

Ra 

Rayleigh  number  Ra  =  g(SNTL3lav 

electrical  conductivity  (fi  1  rrT1) 

Re 

Reynolds  number  Re  =  vfH/rf 

a  thermoelectric  generator,  much  work  has  shown  that  the  heat 
transfer  irreversibilities  between  the  device  and  its  external 
reservoirs  affect  the  performances  of  thermodynamic  processes 
strongly.  Gordon  [32],  Wu  [33],  Agrawal  and  Menon  [34],  Chen  et  al. 
[35,36]  and  Nuwayhid  et  al.  [37,38]  analyzed  the  effect  of  finite-rate 
heat  transfer  between  the  thermoelectric  device  and  its  external 
heat  reservoirs  on  the  performance  of  single-element  thermo¬ 
electric  generators. 

In  practice,  a  commercial  thermoelectric  generator  is  a  multiel¬ 
ement  device,  which  is  composed  of  many  fundamental  thermo¬ 
electric  elements.  Many  researchers  investigated  the  characteristics 
of  multielement  thermoelectric  generators  with  the  irreversibility  of 
finite-rate  heat  transfer,  Joulean  heat  inside  the  thermoelectric 
device,  and  the  heat  leak  through  the  thermoelectric  element.  Esarte 
et  al.  [39]  analyzed  the  influence  of  fluid  flow  rate,  heat  exchanger 
geometry,  fluid  properties  and  inlet  temperatures  on  the  power 
supplied  by  the  thermoelectric  generator.  Chen  et  al.  [40,41] 
analyzed  the  effects  of  external  heat  conductance  and  the  number 
of  elements  on  the  power  and  efficiency  of  a  multielement  gener¬ 
ator.  Yu  and  Zhao  [42]  presented  a  numerical  model  to  predict  the 
performance  of  thermoelectric  generator  with  the  parallel-plate 
heat  exchanger.  Niu  et  al.  [43]  constructed  an  experimental  ther¬ 
moelectric  generator  unit  incorporating  the  commercially  available 


thermoelectric  modules  with  the  parallel-plate  heat  exchanger. 
Behrens  et  al.  [44]  shown  that  the  combustion  of  energy  dense 
liquid  fuels  in  a  catalytic  micro-combustor  can  combust  completely 
for  thermoelectric  applications.  Astrain  et  al.  [45]  studied  the 
influence  of  heat  exchangers’  thermal  resistances  on  a  thermoelec¬ 
tric  generation  system  without  the  internal  factors  of  the  system. 
Hsiao  et  al.  [46]  built  a  one-dimensional  thermal  resistance  model  to 
predict  the  behavior  of  thermoelectric  modules  applied  in  exhaust 
and  radiator  of  an  automobile  without  considering  the  Thomson 
effect  and  the  heat  transfer  through  the  air  gap  of  the  module. 

Reviewing  the  former  literatures  concerning  thermoelectric 
generators,  some  features  can  be  concluded  as  follows: 

(1 )  Researches  on  the  inner  structure  of  a  thermoelectric  generator 
are  without  consideration  of  the  effects  of  external  heat 
transfer  irreversibilities. 

(2)  Researches  on  the  effects  of  external  irreversibilities  of  a  ther¬ 
moelectric  generator  are  without  consideration  of  the  Thom¬ 
son  effect,  inner  structure  and  geometric  dimension  of  the 
thermoelectric  generator. 

(3)  Most  of  the  researches  are  without  consideration  of  the  effects 
of  the  conducting  strips,  ceramic  plate  and  the  heat  transfer 
through  the  air  gap  of  the  module. 
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To  sum  up,  there  is  a  lack  of  a  complete  numerical  model  which 
can  be  applied  to  analyze  the  quantities  of  effects  of  inner  and 
external  factors  on  the  performance  of  a  commercial  thermoelec¬ 
tric  generator  with  heat  exchangers,  and  can  be  applied  to  describe, 
calculate  and  predict  the  practical  performance  of  a  commercially 
available  thermoelectric  generator. 

Taking  into  account  inner  and  external  multi-irreversibilities, 
this  paper  aims  to  establish  a  complete  numerical  model  of 
commercial  thermoelectric  generator  with  finned  heat  exchangers 
by  combining  the  thermodynamics  with  heat  transfer  theory.  Based 
on  the  numerical  model,  the  characteristic  of  a  typical  commercial 
thermoelectric  generator  is  simulated.  The  effects  of  external  irre¬ 
versibilities  on  the  performance  of  the  thermoelectric  generator  are 
analyzed  and  compared  with  the  exo-reversible  model.  The 
numerical  model  and  calculation  method  may  be  applied  to  the 
calculation,  prediction  and  further  optimization  study  for  ther¬ 
moelectric  generators.  The  simulation  results  can  be  used  as 
feasibility  and  effectiveness  reference  by  employing  low-grade 
energy  or  waste  heat  for  power  generation. 


2.  Analysis  of  a  thermoelectric  generator  and  a  complete 
numerical  model 

2.J.  Structure  of  a  thermoelectric  generator  module 

A  typical  commercial  thermoelectric  generator  with  finned  heat 
exchangers  is  shown  in  Fig.  1.  A  thermoelectric  generator  consists  of 
positive-type  (P-type)  and  negative-type  (N-type)  semiconductor 
legs.  The  thermoelectric  semiconductor  legs  are  connected  by  copper 
conducting  strips  in  serials.  A  P-type  and  an  N-type  semiconductor 
leg  compose  a  basic  thermoelectric  element.  The  structure  of  the 
thermoelectric  element  and  the  conducting  strip  is  shown  in  Fig.  2. 
The  length  and  the  side  length  of  the  cross-section  of  a  semi¬ 
conductor  leg  are  respectively  L,  Lb  and  Lb.  La  =  Lb  holds  for  a  typical 
commercial  thermoelectric  module.  The  length  and  the  thickness  of 
the  conducting  strip  are  Lcu  and  ocu,  respectively.  The  number  of 
thermoelectric  elements  is  N.  The  conducting  strips  of  the  thermo¬ 
electric  elements  are  fixed  at  a  thermal  conducting  and  electrical 
insulting  ceramic  plate.  For  a  typical  commercial  thermoelectric 
module,  the  ceramic  plate  is  square.  The  area  and  thickness  of  the 
ceramic  plate  are  respectively  Acp  and  <5cp.  The  module  is  packaging  by 
thermal  insulation  epoxy  resin.  The  thermoelectric  elements  in  most 
commercial  thermoelectric  generator  modules  are  not  closely 
arranged,  that  is,  air  gap  exists  in  the  module.  The  heat  convection  (or 
conduction)  and  radiative  heat  transfer  occur  between  the  two 
ceramic  plates  in  the  area  where  the  thermoelectric  elements  don’t 
occupy.  It  is  assumed  that  the  total  area  of  the  air  gap  is  Ag. 


Heat  exchanger  - 
Ceramic  plate  — 
Thermoelectric  element  - 
Thermal  insulation  - 
Copper  conducting  - 

Heat  exchanger  - 


4  4 


-N- 

Hot  junction 

Cold  junction 

Cold  surface 

4  <0 

r  <0 

-  ^ 

Heat  sink 

-U — [ X2 - 

tricai  current 1 - - — 1  Load  resistance 

Fig.  1.  Structure  of  a  thermoelectric  generator  with  finned  heat  exchangers. 


2.2.  Structure  of  a  finned  heat  exchanger 

A  finned  heat  exchanger  used  for  fluid  heating  or  cooling  is 
shown  in  Fig.  3.  The  base  area,  thickness  and  the  side  length  of 
the  exchanger  are  Ab,  <5exb  and  Lex,  respectively.  The  height,  length 
and  thickness  of  the  fin  are  Hf,  Lf  and  8f,  respectively.  The  width 
of  the  channel  is  r5c.  Lex  —  Lf  holds  for  a  heat  exchanger  with 
square  base. 


2.3.  Thermal  resistance  network 

The  thermal  resistance  network  of  the  whole  device  is  shown 
in  Fig.  4.  All  the  heat  transfers  are  finite-rate  irreversible  heat 
transfer. 

The  hot  and  cold  junction  temperatures  of  the  thermoelectric 
elements  are  Th  and  Tc,  respectively.  The  ceramic  plate  inner  surface 
temperatures  are  Tcpbh  and  Tcpbc.  The  ceramic  plate  exterior  surface 
temperatures  are  Iq^h  and  Tq,jC.  The  fin  base  temperatures  are  Texbh 
and  Texbc-  The  heat  source  (hot  fluid)  and  heat  sink  (cold  fluid) 
temperatures  are  Th  and  Tl,  respectively. 

The  heat  conduction  thermal  resistance  of  a  thermoelectric 
element  is  Rte  =  1  IK  where  K  is  the  thermal  conductance  of 
a  thermoelectric  element.  The  conducting  strip  heat  conduction 
thermal  resistance  are  Rcutc  and  RCuth-  The  heat  convection  (or 
conduction)  and  radiative  thermal  resistance  of  the  air  gap  are  RCVg 
and  Rrdg,  respectively.  The  ceramic  plate  base  heat  conduction 
thermal  resistance  are  Rcph  and  Rcpc.  The  heat  exchanger  base  heat 
conduction  thermal  resistance  are  Rexbh  and  Rexbc.  The  heat 
convection  thermal  resistance  between  the  heat  exchangers  and 
the  heat  reservoirs  are  RCVh  and  Rcvc. 


Fig.  3.  Schematic  diagram  of  a  finned  heat  exchanger. 


3516 


F.  Meng  et  at  /  Energy  36  (201J)  3513-3522 


The  heat  flow  rate  absorbed  from  the  heat  source  to  the  module 
is  Qh-  The  heat  flow  rate  dissipated  from  the  module  to  the  heat 
sink  is  Ql-  The  heat  flow  rate  through  the  hot  and  cold  junctions  of 
the  thermoelectric  elements  are  Qh  and  Qc.  The  heat  flow  rate  of  the 
air  gap  leakage  is  Qg. 

2.4.  Electrical-resistance  network 

The  electrical-resistance  network  of  the  whole  device  is  shown 
in  Fig.  5.  The  electrical  resistance  of  the  P-  and  N-type  semi¬ 
conductor  leg  are  Rp  and  Rn.  The  electrical  resistance  of  the 
conducting  strip  are  RCuh  and  Rcuc.  The  voltage  output,  electrical 
current  output  and  load  resistance  are  U,  1  and  Ri,  respectively. 


2.5.  Energy  equations 

For  a  well-designed  thermal  insulation  packaging  module,  the 
heat  leakage  through  the  surround  of  the  module  can  be  neglected. 
At  steady  work  state,  the  temperature  distribution  of  the  air  gap 
is  the  same  with  the  thermoelectric  elements,  so  the  heat  transfer  of 
the  whole  device  can  be  treated  as  one-dimensional  heat  transfer. 

The  infinitesimal  as  shown  in  Fig.  6  is  considered.  According 
to  the  non-equilibrium  thermodynamics,  the  inner  effects  of  the 
thermoelectric  elements  include  Seebeck  effect,  Fourier  effect, 
Joule  effect  and  Thomson  effect.  The  increment  rate  of  inner  energy 
of  the  infinitesimal  is  zero  at  steady-state,  so  one  can  obtain  the 


Fig.  5.  Schematic  diagram  of  the  electrical  resistance  network. 


energy  conservation  equation  of  P-type  semiconductor  leg  as 
follows  [1,2] 

OKin  -  Qltout  +  Q|  +  Qp  =  0  (1 ) 

where  Ogm,  Qicout,  Q|  and  Qp.  are  the  Fourier  heat  input,  Fourier  heat 
output,  generated  Joule  heat  and  Thomson  heat,  respectively.  Eq. 
(1 )  can  be  written  as 

A(Tp  +  dTp)  Mp -fMP+  (e-/P)2^f + VdTP  =  0  (2) 

where  kp,  ap,  pp,  Ap,  Lp,  Tp,  ejp  are  respectively  the  thermal 
conductivity,  electrical  conductivity,  Thomson  coefficient,  cross- 
section  area,  length,  temperature  and  electrical  current  density  of 
the  P-type  semiconductor  leg.  Reforming  Eq.  (2)  and  making  same 
analysis  on  the  N-type  semiconductor  leg,  one  can  obtain  the  heat 
conduction  differential  equations  (3)  and  (4),  with  boundary 
conditions  (5)  and  (6)  as  follows  [11] 


&l  e- 

+ 

+ 

O 

(3) 

(4) 

Tp(0)  =  Tn(0)  =  Tc 

(5) 

Tp(Lp)  =  Tn(Ln)  =  Th 

(6) 

where  kn,  on,  fin,  An,  Tn,  ejn  respectively  are  the  thermal  conduc¬ 
tivity,  electrical  conductivity,  Thomson  coefficient,  temperature 
and  electrical  current  density  of  the  N-type  semiconductor  leg. 

Taking  into  account  the  effect  of  temperature  dependence  of 
thermoelectric  properties,  kp,  op  and  pp  are  function  of  Tp;  while  kn, 
on  and  fin  ate  function  of  Tn.  However,  such  a  differential  equation 
cannot  be  solved  analytically.  Replacing  k,  a  and  q  with  the  mean 
value  k,  o  and  p  approximately  gives  approximation  of  Eqs.  (3)  and 
(4)  as  follows 

(7) 

r  .  d T2  _  dTn  / 2 

KiAi-r4  -  Pnf—j —  +  -  ,  =  0 

dx2  dx  On  An 

(8) 

where  kp  =  kp|r=(rh+iy/2.  °p  =  ap\T=(Th+Tc)/2’  Mp  =  MpIt=ct]1+tc)/2 
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for  P-type  semiconductor  leg  and  kn  —  kn\T=^+Tcy2< 
On  =  on\T=(Th+Tc)/2'  Pn  =  MnlT=(Th+rc)/2  f°r  N-type.  Alteration 
method  along  with  the  following  fitted  formulas  Eqs.  (53)-(55) 
are  adopted  to  determine  the  junction  temperatures  Th  and  Tc. 

The  total  heat  flow  rates  through  the  hot  and  cold  junctions  are 

«*  -  N((«ph  -  «■*)  V  +  +  W£|„J  (9) 

a  -  N((«p,-«„)rci+Mp^ja+Vi«^|j=0)  (10) 

where  ap  and  an  are  the  Seebeck  coefficient  of  the  P-  and  N-type 
semiconductor  leg,  and  the  subscript  h  and  c  represent  the  hot  and 
cold  side,  respectively. 

The  heat  flow  rates  through  the  conducting  strips  are  given  by 

Qhj  s  Jv(rcpbh  -  Th) / Rcmh  (H) 

Qci  =  N(TC  -  Tcpbc) /f W  (12) 

It  is  assumed  that  the  total  thermal  resistance  between  the  hot 
side  ceramic  plate  inner  surface  and  the  heat  source  is 
Reh  =  Rcvh  +  Pexbh  +  Pcph  (as  shown  in  Fig.  4);  the  total  thermal 
resistance  between  the  cold  side  ceramic  plate  inner  surface  and 


the  heat  sink  is  Rec  =  RCvc  +  Rexbc  +  RCpc-  Then,  one  has 

Qh  =  (TH-rcpbh)/Reh  (13) 

Ql  =  (j^cpbc  Ti)  /Re  c  (14) 

The  heat  flow  rate  through  the  air  gap  (gray  arrows  in  Fig.  1 )  Qg 
can  be  expressed  as 

Qg  =  Qcvg  +  Qrdg  (15) 

Qcvg  =  (Vcpbh  -  Icpbc)  / Revs  (15) 

Qrdg  =  (Tc4pbh  I'cpbc)  /  Rrdg  (12) 


where  Qcvg.  Qrdg.  Rcvg  and  Rrdg  represent  the  convection  heat  rate, 
radiative  heat  rate,  heat  convection  resistance  and  heat  radiative 
resistance,  respectively. 

The  relations  of  heat  flow  rates  can  be  expressed  as  follows: 


Qh  —  Qh  +  Qg 

(18) 

Qh  =  Qhi 

(19) 

Qc  =  Qd 

(20) 

Ql  =  Qc  +  Qg 

(21) 

3.  Solution  of  the  model 

3.1.  Solution  of  the  thermoelectric  element 

The  temperature  distribution  of  the  P-  and  N-type  semi¬ 
conductor  leg  can  be  obtained  by  solving  Eqs.  (7)  and  (8)  as  follows: 


7p(x)  =  Tc  -  Fpx  +  Th(eJX/_Ff(e^  ~  l)  (22) 


MX)  =  TC  +  Fnx  +  t23) 


where  wp  =  pp//(kpAp),  Fp  =  I/(ap/cpA),  <o„  =  nnI/(knAn)  and 
Fn  =  I/QfnhA). 

For  maximum  figure  of  merit  of  the  thermoelectric  element 
Z  =  u2I(KR),  the  geometry  dimension  of  the  thermoelectric  element 
and  the  physical  property  of  the  material  should  satisfy  the 
following  equation  [1,2] 


Afc_knOn 
Aiq  kp(T 


To  reduce  the  cost  of  manufacture,  the  P-  and  N-type  semi¬ 
conductor  legs  are  made  with  same  geometry  dimensions  i.e. 
Ap  =  An  =  A  and  Lp  =  Ln  =  L.  So  one  has  kD<rJ(kpap)  -  1  by  Eq.  (24). 
Similar  doped  alloys  are  adopted  to  make  P-  and  N-type  semi¬ 
conductor  legs.  That  is  <rp  =  <jn  =  a,  kp  =  kn  =  k,  ap  =  -a„,  /ip  =  -fin. 
According  to  Taylor’s  formula,  when  |x|  <s  1,  e*  =  1  +  x  holds  true. 
Based  on  above  assumption,  Eqs.  (9)  and  (10)  can  be  approximated 
as  follows  at  the  case  of  fil/K  1 


Qh  =  N[abrrh  +  K(Th  -  Tc)  -  0.5 12R  -  0.5 ]U(Tb  -  Tc)]  (25) 


Qc  =  N [ac/Tc  +  K(Tb  -  Tc)  +  0.5 12R  +  0.5 JlI(Tb  -  Tc)]  (26) 

where  ab  =  ahp  -  ahn  and  ac  =  acp  -  acn  are  the  Seebeck  coefficient 
of  the  thermoelectric  elements  at  hot  and  cold  side. 
Jl  =  ftp  -  Jin  =  2pp,  K  and  R  are  respectively  the  total  Thomson 
coefficient,  thermal  conductance  and  electrical  resistance  of 
a  thermoelectric  element.  alT,  KAT,  I2R  and  filAT  are  the  rates  of 
Peltier  heat,  Fourier  heat,  Joule  heat  and  Thomson  heat,  respec¬ 
tively.  K  and  R  are  given  by 

K  =  Kp  +  Kn  =  kpAp/Lp  +  knAn/Ln  =  2  kA/L  (27) 


R  =  Rp  +  Rn  =  Lp/(opAp)  +  Ln/(onAn)  =  2  L/(oA)  (28) 

The  power  output  is  the  difference  between  the  absorbed  and 
dissipated  heat  flow  rate  and  results  in 

P  =  Qh  -  Qc  =  NI[ah  Tb  —  ac  Tc  —  IR  -  p(Th  -  Tc)]  (29) 

Taking  the  Joule  heat  generated  by  the  conducting  strips  as  the 
thermoelectric  elements  generated  and  replacing  R  with 
Rt  =  R  +  Rcu  gives  the  modified  power  output  as  follows 


P  =  NI[ah  Th  -  ac  Tc  -  IRt  -  p(Th  -  Tc)] 

The  thermal  efficiency  r]  =  P/Qh  is  given  by 
M[ah  Th  -  ac  Tc  -  IRt  -  p(Th  -  Tc)] 

(Th  -  Tcpbh)  / Reh 

Taking  the  heat  source,  heat  sink,  thermoelectric  generator  and 
heat  exchangers  as  an  isolated  system,  the  entropy  production 
rate  of  the  system  is  the  entropy  increase  of  the  system,  and 
results  in 

Sg  =  ASh  +  ASl  +  ASteg  =  -Qh/Th  +  Ql/Tl  (32) 


(30) 

(31) 
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where  ASh,  ASl  and  ASteg  are  the  entropy  increase  of  the  heat 

source,  heat  sink  and  the  thermoelectric  generator,  respectively.  ^  =  2A/V/Acp  (44) 

ASteg  =  0  because  the  thermoelectric  generator  is  in  steady.  The  electrical  resistance  of  a  conducting  strip  is  given  as  follows 


3.2.  Calculation  method  of  the  thermal  and  electrical  resistances 

The  heat  conduction  thermal  resistance  of  a  conducting  strip 
and  a  ceramic  plate  are  as  follows: 


RcP  =  5cp/(kcpAcp)  (34) 

The  heat  conduction  thermal  resistance  of  the  heat  exchanger 

Kexb  =  <W(Wlb)  (35) 

According  to  the  heat  transfer  theory  [47^19],  the  heat 
convection  thermal  resistance  between  the  finned  heat  exchanger 
and  the  heat  reservoir  is 


Re 


hcjAeff 


(36) 


where  ho,  and  Aefr  are  respectively  the  coefficient  of  heat  convec¬ 
tion  and  the  effective  heat  transfer  area,  hcv  can  be  calculated  by 


hcv 


Nu  kf 

~W 


(37) 


where  Nu,  kf  and  Hf  are  the  Nusselt  number  of  the  flow,  thermal 
conductivity  of  the  fluid  and  the  height  of  the  fins,  respectively.  Nu 
can  be  calculated  by  [49] 


Nu  =  0.664Re1/2Pr1/3,  Re  <  5  x  105  (Laminar  flow) 

Nu  =  0.037Re4/5Pr1/3,  Re  >  5  x  105  (Turbulent  flow) 

where  Prf  and  Re  are  the  Prandtl  number  of  the  fluid  and  the 
Reynolds  number  of  the  flow,  respectively.  They  can  be  calculated 
by 


Ptf 


bf  =  cpt<Pf 
a{  kf 


(39) 


Re  =  VfHf/Vf  (40) 

where  vf,  at,  cPf,  kf,  cpt,  pt  and  Vf  are  the  Kinematic  viscosity,  thermal 
diffusivity,  constant  pressure  heat  capacity  ratio,  thermal  conduc¬ 
tivity,  kinetic  velocity,  density  and  the  viscosity  of  the  fluid, 
respectively.  The  effective  heat  transfer  area  is  given  by  [47] 


Aeff  =  Vf  (2NfLfHf)  +  (Nf  -  l)  8cLf  (41 ) 

where  the  fin  efficiency  Vf  is  given  by 


To  describe  the  geometry  of  the  thermoelectric  module  air  gap, 
a  ratio  is  introduced  from  Ref.  [13]  and  defined  as  module  packed 
density  6 


Rcu  =  JL  =  1  ,?&  ■  ia(i  -  B)/e  =  l  (45) 

ffcu  ^cu^b  ^CU  iHilji,  ^cu  Ifeti. 

where  a€U  and  5CU  are  the  electrical  conductivity  and  thickness  of 
the  conducting  strip,  respectively. 

The  air  gap  of  the  module  can  be  considered  as  a  restricted 
horizontal  interlayer.  There  are  two  cases:  when  the  underside  is 
cold  side,  the  heat  transfer  is  heat  conduction;  when  the  underside 
is  hot  side,  the  heat  transfer  mode  is  dependent  on  the  Ra  number 
(the  product  of  Cr  and  Pr)  which  is  given  by  [47] 


Ra 


g/3akATL3  uair 

uair  aair 


gfeirATI3 

aajrl’air 


(46) 


where  /3air.  aair,  fair  and  g  are  the  coefficient  of  cubical  expansion, 
thermal  diffusivity  and  kinematic  viscosity  of  air,  and  gravitational 
acceleration,  respectively. 


(1)  When  Ra  <  1700,  the  heat  transfer  is  heat  conduction  and  then 
the  thermal  resistance  is 


^edg 


L  L 

-  ACpfcair(l  -  6) 


(47) 


(2)  When  Ra  >  1700,  the  heat  conduction  develops  to  natural 
convection  and  then  the  thermal  resistance  is 


o  _  1 

Kcvs  -  u — 2" 

ricvg/ig 

where 

/lev g  Nu  kjif/L 


Nu  =  0.59  Ra04, 
Nu  =  0.212  Ra1/4, 
Nu  =  0.061  Ra1/3, 


1.7  x  103  <  Ra  <  7.0  x  103 
7.0  x  103  <  Ra  <  3.2  x  105 
Ra  >  3.2  x  105 


(48) 


(49) 

(50) 


The  radiative  heat  flow  rate  between  two  parallel  plates  with 
same  area  is  given  by 


Qrd 


^(t4-T4)  ta_ta 

rdl/ei  +  l/e2  -  1  Rrdg 


(51) 


where  ei,  e2,  Ti,  T2  and  Ar d  are  the  blackness,  temperature  and  area 
of  the  two  plates,  respectively,  at,  is  the  Stefan-Boltzmann 
constant.  So  one  can  derive  the  radiative  thermal  resistance  of  the 
air  gap  as  follows 


Rrdg  — 


l/ecp  +  l/£cp  —  1  _  2/ eCp  —  1 

abAs  ~  abAcp(l  -  6) 


where  eCp  is  the  blackness  of  the  ceramic  plates. 


(52) 


Table  1 

Geometric  dimensions  of  the  thermoelectric  module  [51]. 

A  (mm2)  I  (mm)  <5CU  (mm)  <5cp  (mm)  Acp  (mm2) 

~  L9  02  09  29.7  x  29.7  ^ 
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Table  2 

Physical  properties  of  conducting  strip,  ceramic  plate,  and  heat  exchangers  [52], 
Pea  (Qm)  tcutWm-’r1)  kcpfWm-'r1)  ecp  (W  m-1  IT1) 

1.7  x  IQ-8  386  353  09  204 


4.  Simulation  and  comparison 

The  characteristic  of  the  thermoelectric  generator,  including 
voltage  output,  power  output,  efficiency,  entropy  production  rate, 
temperature  and  heat  flow  rate  features  is  investigated  by 
numerical  simulation.  Water  represents  a  growing  concern  for 
meeting  future  power  generation,  especially  thermoelectric 
power  generation  needs  [50],  Almost  all  of  the  energy  can  be 
absorbed  by  the  water  to  convert  to  available  energy.  So  hot  water 
at  60—100  °C  and  cold  water  at  27  °C  (room  temperature)  are 
employed  as  heat  source  and  sink  of  the  generator.  A  typical 
commercial  thermoelectric  generator  made  by  Ferrotec  (China) 
[51]  (8001/127/040B)  which  consists  of  127  thermoelectric 
elements  is  chosen  for  this  simulation.  The  simulation  results  can 
be  used  for  reference  by  employing  low-grade  energy  or  waste 
heat  for  power  generation. 

The  physical  properties  of  the  commercially  available  material 
by  Melcor  [52]  used  for  this  simulation  are  shown  as  follows 

ap  s  |22224.0+  930.6T-  0.9905T2)lCT9  VK1  (53) 

p  =  (5112.0+  163.4T  +  0.6279T2)lCT10  Qm  (54) 

k=  (62605.0  -  277.7T  +  0.4131T2)  10  4  Wm1  KT1  (55) 

where  ap,  p  and  k  are  the  Seebeck  coefficient,  electrical  resistivity 
and  thermal  conductivity,  respectively.  T  =  (Tj,  +  Tc)/2  is  the  mean 
temperature  of  the  hot  and  cold  junction. 

The  geometric  dimensions  of  the  thermoelectric  module  are 
shown  in  Table  1  [51],  The  physical  properties  of  the  conducting 
strip,  ceramic  plate  and  heat  exchanger  are  shown  in  Table  2  [52], 
The  geometric  dimensions  of  heat  exchangers  are  shown  in  Table  3 

[47] ,  The  physical  properties  and  flow  parameters  of  air  and  water 

[48]  are  shown  in  Tables  4  and  5,  respectively. 

Fig.  7  shows  the  temperature  difference  between  the  hot  and 
cold  junction  AThC  =  Th  —  Tc  versus  electrical  current  output.  It  can 
be  seen  that  the  temperature  difference  is  linear  and  decreasing 
function  of  the  electrical  current  output.  Actually,  the  linear  rela¬ 
tionship  is  approximate  because  there  are  factors  (such  as  heat 
convections  and  radiative  heat  transfer)  which  caused  non¬ 
linearity.  However,  they  affected  the  linearity  little.  The  tempera¬ 
ture  difference  is  affected  by  the  electrical  current  output  little,  but 
affected  by  the  hot  water  temperature  strongly. 

Fig.  8  shows  the  voltage  output  versus  electrical  current  output. 
The  curves  plotted  by  dotted  lines  in  this  figure  (and  following 
Figures)  are  calculated  based  on  the  exo-reversible  model,  i.e.,  all  the 
external  irreversibilities  are  neglected.  It  can  be  seen  that  the  output 
voltages  are  linear  and  decreasing  functions  of  the  electrical  current 
output.  The  slopes  of  the  curves  are  the  same.  That  is  for  any  output 


Table  3 

Geometric  dimensions  of  heat  exchangers  [47], 

Ab  (mm2)  6b  (mm)  5f  (mm)  6C  (mm)  Hf  (mm)  Nf 

29.7  x  29.7  4  Z7  2J  20  6~” 


Table  4 

Parameters  of  air  [48], 

(^(Wm-’K-1)  Prair  feirOt-1)  rpn,  (kgm-1  s'1)  i'air(m2s“1) 

2.57  x  IQ-2  0.713  3.43  x  10~3  1.82  x  10~5  1.52  x  10~5 


current,  the  external  irreversibilities  decrease  the  output  voltage 
with  a  same  value  10%  and  the  heat  reservoir  temperature  difference 
does  not  affect  the  slopes  of  the  curves.  Specially,  when  the  hot  water 
temperature  is  100  °C,  the  open  circuit  voltage  Uopen  =  1.80  V  (/  =  0) 
and  the  short  circuit  current  fSh0rt  =  0.28  A(U=  0). 

Fig.  9(a)  and  (b)  shows  the  power  output  and  efficiency  versus 
electrical  current  output,  respectively.  It  can  be  seen  that  P—I  and 
t)—I  curves,  in  Fig.  9(a)  and  (b),  are  both  parabolic-like,  which  means 
there  is  one  optimal  current  output  f  corresponding  to  the 
maximum  power  output  Pmax  and  another  optimal  current  output 
f1  corresponding  to  the  maximum  efficiency  r/max.  The  external 
irreversibilities  decrease  power  output  and  efficiency  by  17%  and 
5%,  respectively.  When  the  hot  water  temperature  increases,  the 
maximum  power  and  maximum  efficiency  along  with  the  corre¬ 
sponding  optimal  currents  increase.  When  the  hot  water  temper¬ 
ature  is  100  °C,  the  maximum  power  output  of  0.13  W  and  the 
maximum  efficiency  of  0.87%  are  available  from  the  generator. 

Fig.  10  shows  the  entropy  production  rate  versus  electrical 
current  output.  It  can  be  seen  that  the  entropy  production  rate  is 
increasing  function  of  the  electrical  current  output.  When  J  — >  0, 
the  entropy  production  rate  tends  to  the  minimum.  The  reason  is 
that  the  entropy  production  rate  by  the  Joule  heat  tends  to  zero,  all 
the  entropy  production  is  generated  by  inner  and  external  heat 
transfer  which  can  be  expressed  by 
Sgj_o  -  Qh(1/Tl  -  1/Th)  =  Ql(1/Tl-  1/Th).  It  should  be  noted 
that  the  external  irreversibilities  decrease  the  entropy  production 
rate  for  fixed  electrical  current  output.  On  the  one  hand,  the 
external  irreversibilities  cause  external  entropy  production  but  the 
lower  temperature  difference  between  the  hot  and  cold  junction  of 
the  thermoelectric  elements  decreases  the  inner  entropy  produc¬ 
tion  rate.  On  the  other  hand,  the  external  irreversibilities  decrease 
the  rate  of  heat  transfer  from  the  heat  source  to  the  heat  sink.  When 
the  external  thermal  resistance  tends  to  infinite,  the  rate  of  heat 
transfer  and  then  the  entropy  production  tends  to  zero.  Certainly, 
the  decrease  of  entropy  production  rate  is  at  a  cost  of  the  decrease 
of  power  output  and  efficiency.  On  the  whole,  the  entropy  gener¬ 
ation  minimization  is  not  suitable  for  the  electrical  current  opti¬ 
mization  of  a  thermoelectric  generator.  Only  power  and  efficiency 
maximization  are  useful  and  possible  in  this  case. 

Fig.  11(a)  and  (b)  shows  the  surface  temperatures  versus  elec¬ 
trical  current  output.  It  can  be  seen  that  all  the  hot  side  tempera¬ 
tures  decrease  while  all  the  cold  side  temperatures  increase  with 
the  increase  in  the  electrical  current  output.  So  the  temperature 
difference  between  the  hot  and  cold  junctions  decreases  (Fig.  7). 
Comparing  the  temperature  difference  caused  by  different  parts, 
the  following  equation  holds:  AT„  >  ATcp  >  ATexb  >  ATCU,  where 
ATCV,  ATcp,  ATeXb  and  ATCU  are  the  temperature  differences  caused  by 
the  fins  heat  convection,  ceramic  plate,  heat  exchanger  base  and 
conducting  strips,  respectively.  The  temperature  differences  reflect 
the  different  effect  of  the  thermal  resistance  on  the  performance  of 
the  device. 


Table  5 

Parameters  of  water  [48], 

U,  (W  m-1  K-1)  Prwater  v  (m  s-1)  TH("C)  TL(“C)  AT(°C) 

5.99  x  10-2  Z02  05  60-100  27  33-73 
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Fig.  7.  Temperature  difference  between  the  hot  and  cold  junctions  versus  electrical  t) 

current  output.  TH  and  7j.  are  the  temperatures  of  the  hot  and  cold  water,  respectively. 


Fig.  12(a)  and  (b)  shows  the  heat  flow  rates  versus  electrical 
current  output.  It  can  be  seen  that  the  heat  flow  rates  can  be  ^ 
classified  into  two  groups  according  to  the  quantity  of  the  heat  flow 
rates.  Heat  flow  rates  in  one  group  are  of  large  quantity  (Fig.  12(a))  o' 

and  others  in  another  group  are  of  small  quantity  (Fig.  12(b)).  Qh,  Ql.  .1 

Qh.  Qc>  Qj,  Qa  and  are  increasing  functions  of  current  output.  The 
inner  heat  leakage  Qk  decreases  because  of  the  decrease  of  A7hc  W 
(see  Fig.  7).  It  should  be  noted  that  the  Thomson  heat  is 
considerable  and  cannot  be  neglected.  Among  the  heat  flow  rates 
caused  by  inner  effects  i.e.  the  rates  of  Peltier  heat,  Fourier  heat, 

Joule  heat  and  Thomson  heat,  the  Fourier  heat  leakage  is  the  main 
loss.  The  inner  heat  leakage  doesn’t  affect  the  power  output  but 
increase  the  absorbed  heat  and  then  decrease  the  efficiency.  The 
Seebeck  power  is  a  linear  function  of  current  output  approximately, 
but  the  Joule  heat  is  a  quadratic  function  of  electrical  current 
output.  This  is  the  reason  why  there  is  an  optimal  current  corre¬ 
sponding  to  the  maximum  power  output.  The  calculation  shows 
that  the  heat  transfer  through  the  air  gap  is  heat  conduction  and 
radiative  heat  transfer  because  the  heat  transfer  cannot  develop  to 
natural  convection  heat  transfer.  The  narrow  space  and  low 
temperature  difference  are  the  main  factors.  However,  the  heat 


Electrical  current,  1(A) 


Fig.  9.  (a)  Power  output  versus  electrical  current  output.  TH  and  TL  are  the  tempera¬ 
tures  of  the  hot  and  cold  water,  respectively,  (b)  Efficiency  versus  electrical  current 
output.  Th  and  TL  are  the  temperatures  of  the  hot  and  cold  water. 


Electrical  current,  1(A) 


Electrical  current,  1(A) 


Fig.  8.  Electrical  voltage  output  versus  electrical  current  output.  TH  and  TL  are  the  Fig.  10.  Entropy  production  rate  versus  electrical  current  output.  TH  and  Tl  are  the 
temperatures  of  the  hot  and  cold  water,  respectively.  temperatures  of  the  hot  and  cold  water,  respectively. 
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Fig.  11.  (a)  Surface  temperatures  versus  electrical  current  output.  Th,  Tcpbh,  Tcpib,  Texbh 
are  thermoelectric  element  junction  temperature,  ceramic  plate  inner  and  exterior 
surface  temperature,  and  fins  base  temperature  of  the  hot  side,  respectively.  TH  and  TL 
are  the  temperatures  of  the  hot  and  cold  water,  respectively,  (b)  Surface  temperatures 
versus  electrical  current  output.  Tc,  Tc^c  T^,  Texbc  are  thermoelectric  element  junc¬ 
tion  temperature,  ceramic  plate  inner  and  exterior  surface  temperature,  and  fins  base 
temperature  of  the  cold  side,  respectively.  TH  and  TL  are  the  temperatures  of  the  hot 


flow  rate  through  the  air  gap  Qg  is  considerable  and  even  larger 
than  the  power  output  generally.  The  heat  flow  rate  through  the  air 
gap  decreases  with  the  increase  in  electrical  current  output  because 
of  the  decrease  of  AThc  (Fig.  7). 


-  Irreversible 

- Exo-reversible 


Electrical  current,  1(A) 


Fig.  12.  (a)  Heat  flow  rates  versus  electrical  current  output.  Qh  is  the  total  absorbed 
heat  rate  input  to  the  thermoelectric  generator,  Ql  is  the  total  rejected  heat  rate  output 
to  the  thermoelectric  generator,  Qh  is  the  absorbed  heat  rate  input  to  the  thermo¬ 
electric  elements,  Qc  is  the  rejected  heat  rate  output  to  the  thermoelectric  elements, 
Qk  =  NK(Th  -  Tc)  is  the  inner  heat  leakage,  and  P  =  QH-Qi.  =  Qh-Qcis  power  output. 
Th  and  TL  are  the  temperatures  of  the  hot  and  cold  water,  (b)  Heat  flow  rates  versus 
electrical  current  output.  Q,  =  NaI(Th  -  Tc )  is  the  Seebeck  power  rate,  Qj  =  NI2R  is  the 
Joule  heat  rate,  Qg  =  JVpf(Tb  -  Tc)  is  the  Thomson  heat  rate,  and  P  =  Q;,  Qj  Q;1  is  power 
output,  respectively.  Qg  =  Qnvg  +  Qrdg  is  the  heat  flow  rate  through  the  air  gap,  where 
Qcvg  and  Qrdg  represent  the  heat  conduction  rate  and  heat  radiative  rate,  respectively. 
Th  and  Tl  are  the  temperatures  of  the  hot  and  cold  water,  respectively. 


5.  Conclusions 

A  complete  numerical  model  with  inner  and  external  multi¬ 
irreversibilities  of  commercial  thermoelectric  generator  with  fin¬ 
ned  heat  exchangers  is  established,  in  which  physical  properties, 
geometric  dimensions,  and  flow  parameters  are  all  considered. 

Applying  the  numerical  model  to  a  typical  commercial  ther¬ 
moelectric  generator,  the  characteristic  of  the  thermoelectric 
generator  is  simulated.  Hot  water  at  60—100  °C  and  cold  water  at 
27  °C  (room  temperature)  are  employed  as  heat  source  and  sink  of 
the  generator.  It  is  found  that 

(1)  The  temperature  difference  between  the  junctions  and  voltage 
output  is  a  linear  function  of  the  electrical  current  output 
approximately.  P—  I  and  77 — /  curves  are  both  parabolic-like. 
When  the  hot  water  temperature  is  100  °C,  the  open  circuit 


voltage  is  1.80  V,  the  short  circuit  current  is  0.28  A,  the 
maximum  power  output  is  0.13  W  and  the  maximum  efficiency 
is  0.87%,  respectively.  The  entropy  generation  minimization  is 
not  suitable  for  the  electrical  current  optimization  of  a  ther¬ 
moelectric  generator. 

(2)  The  Fourier  heat  leakage  is  the  main  loss  among  the  losses 
caused  by  inner  effects.  The  heat  transfer  through  the  air  gap  is 
heat  conduction  and  radiative  heat  transfer.  The  heat  flow 
through  the  air  gap  and  Thomson  heat  is  considerable  and  even 
larger  than  the  power  output  generally.  Sorting  by  the  quantity 
of  effect  on  the  performance  of  the  device,  the  maximum  is  the 
thermal  resistance  of  the  heat  convection  between  the  heat 
exchanger  and  the  fluid,  followed  by  the  heat  leakage  through 
the  air  gap,  ceramic  plates,  and  thermal  resistance  of  the  heat 
exchanger  base,  the  minimums  are  the  thermal  and  electrical 
resistance  of  the  conducting  strips. 


